ORBITALLY STABLE STANDING WAVES FOR THE 
ASYMPTOTICALLY LINEAR ONE-DIMENSIONAL NLS 
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Abstract. In this article we study the one-dimensional, asymptotically lin- 
ear, non-linear Schrodinger equation (NLS). We show the existence of a global 
smooth curve of standing waves for this problem, and we prove that these 
standing waves arc orbitally stable. As far as we know, this is the first rig- 
orous stability result for the asymptotically linear NLS. We also discuss an 
application of our results to self-focusing waveguides with a saturable refrac- 
tive index. 



1. Introduction 

In this article we study the one-dimensional nonlinear Schrodinger equation 

id^ + d 2 xx ^ + /(x, |V| 2 )-0 = (NLS) 

for ip = ip(t, x) : [0,oo)xR^ C. We suppose that / e C 1 (M x K + , M), and satisfies 
the following assumptions. 

(AO) f(x, 0) = for all and we have 

lim fix, s) = 0, uniformly for x € R, (1-1) 

lim f(x,s) = 0, uniformly for s 0. (1-2) 

\x\— >oo 

(AL) There exists e C(R) n L°°(R) such that 

lim f(x,s) = foo(x), uniformly for x G M. (1-3) 

s— >oo 

Equation (|NLS[) with the assumption (AL) is usually referred to as the asymp- 
totically linear nonlinear Schrodinger equation. The structure of the nonlinearity 
allows one to look for standing wave solutions 

i/j(t, x) = e iXt u(x), where A > and u : R ->• R. 

As we shall see from the stability analysis in Section [3j solutions of this type enjoy 
remarkable properties with respect to the general dynamics of (|NLSj) . Using this 
Ansatz to solve (jNLSp yields the stationary equation 

u" + f(x, u 2 )u = Xu. (SNLS) 

This second order ordinary differential equation will be interpreted as a nonlinear 
eigenvalue problem, which will be addressed via bifurcation theory. 
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By a solution to (jSNLSj) will be meant a couple (A, u) £ M x iJ 1 (R) satisfying 
(jSNLSj) in the sense of distribution. Note that, since f(x, 0) = 0, we have a line of 
trivial solutions, {(A, 0) : A £ 1} C 1 x H 1 (M.). We will prove that there exist non- 
trivial solutions, bifurcating from the line of trivial solutions. Under appropriate 
hypotheses on /, in particular assuming that f(x, s) is positive and even in x, we 
will obtain a smooth curve of positive even solutions of (jSNLSj) . 

{(A, u{\)) : A e (0, Aoo)} Clx H^R). (1.4) 

The number Aoo G (0, oo) will be characterized as the principal eigenvalue of the 
linear problem 

u" + foo(x)u = Xu, (1.5) 

known as the asymptotic linearization. 

Additional properties of the solutions will be proved, in particular their bifurca- 
tion behaviour: 

lim ||u(A)|| ff i (R) = and lim ||u(A)|| H i (K ) = oo. 

The curve (|1.4[) will be obtained, first by a local bifurcation analysis near u = 0, then 
by analytic continuation and an asymptotic analysis as ||u(A)||#i — > oo. Bifurcation 
from the line of trivial solutions is difficult in the present context since, under 
assumption the linearization of (jSNLSj) at it = is 

u" = Xu, (1.6) 

and has purely continuous spectrum. 

In Subsection 12. 1[ prescribing the behaviour of f(x,s) as s — ¥ 0, we will obtain 
bifurcation of positive solutions of (jSNLSj) . from the bottom of the continuous 
spectrum of (jl.6| (i.e. from A = 0), by applying a fairly involved perturbation 
analysis, introduced in 0]. This local result is contained in Theorem 12.31 

Under additional assumptions (in particular assuming that f{x, s) is positive and 
even in i g R), a global analysis is carried out in Subsection 12.21 where we show 
that the local branch given by Theorem 12.31 can be extended in a smooth manner 
to the curve (|1.4j) . The global continuation relies on the non-degeneracy of positive 
even solutions of (jSNLSj) given by Lemma 12.91 The asymptotic behaviour of the 
branch as |ju(A)||#i — > oo (and A — > Aoo) follows by the asymptotic bifurcation 
analysis that was developped in Q using topological arguments. Our global result 
is Theorem 12. 101 

In (where we only considered (jSNLSj) for x > 0) we obtained global bifurcation 
from (Aoo, oo) in R x H 2 (0, oo). More precisely, we obtained a global connected set 
of solutions (A, u) such that ||M|| ff 2 — > oo as A — > Aoo- This was later extended in [l(| 
to a higher dimensional version of (jSNLSj) . under fairly weak hypotheses, providing 
a strong existence result for positive solutions of (jSNLSI) with assumption (AL). 
Previous contributions on this problem — see e.g. ^TEE EH I27W30I and the 
references in these papers — mostly used variational methods to prove existence 
of solutions for (jSNLSj) - (AL). under various assumptions, typically stronger than 
those required by the topological approach. 

However, due to technical restrictions of the method used in 0], we were not 
able to continue the bifurcating branch down to the line of trivial solutions. Sec- 
tion [2] below completes the discussion initiated in by showing that, under more 
restrictive assumptions (in particular, symmetry and monotonicity conditions on 
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f(x, s), as well as a precise asymptotic behaviour as s — > 0), there exists a smooth 
curve of solutions, connecting (0,0) to (Aoo,oo) in R x H 1 (M.). 



A global solution curve was obtained by Jeanjcan and Stuart 15| under similar 
hypotheses. However, they consider an additional, non-trivial, linear potential in 
(ISNLSp . In our notation, this amounts to assuming that fo(x) := /(x, 0) ^ 
instead of f(x, 0) = in (AO). In this case, the linearization at u = has the form 

u" + f a (x)u = Xu. (1.7) 

Under appropriate assumptions on f (e.g. fo is a 'bump'), the linear problem ()1.7[) 
has a principal eigenvalue, from which bifurcation can be obtained via standard bi- 
furcation theory. Global continuation can then be obtained by arguments similar 
to those of Subsection l2.2l The authors of [TBI also discussed the case of an asymp- 
totically linear nonlinearity, as in (| 1 .3[) . The main difference in the present context 
is that bifurcation at u — occurs from the bottom of the continous spectrum of 
the linearization (|1.6j) . 

In Scction[3]we will consider the standing wave solutions ip\(t,x) = e lXt u(X)(x) 
of (|NLS|) corresponding to the solutions (TO} of (jSNLSj) . We will prove that they 
are orbitally stable amongst the set of solutions ■(/>(£, a;) <G C([0,oo),if 1 (R,C)). This 
result, Theorem 13.81 is based on the general theory of orbital stability for Hamil- 
tonian systems, see |l2l[26l|. Given a standing wave ip\ Q , it follows from the theory 
that, under appropriate conditions on the spectrum of the linearization of (jNLSj) 
at V'Ao: this standing wave is orbitally stable if the mapping A — > \\u(X)\\ 2 L2 is in- 
creasing at the point A — Xq. Since we have a smooth curve of solutions, this can 
be obtained by checking that 



_d_ 

dA 



|KA)||£ 2 > 0. (1.8) 

A = An 



Our bifurcation analysis in Section [2] shows that ||zt(A) || x, 2 — > as A — > 0, and so 
(|1.8|) must hold for some Ao £ (0,Aoo). Hence, by continuity, we need only check 
that 

-£-||u(A)||£ a ^0, forall AG (0,Aoo). 

This is done in Subsection 13.21 using an integral identity that was first derived in 
fl(| to study orbital stability along the solution curve obtained in (lEij] . However, 
the authors of [3] were not able to deal with the asymptotically linear case. In 
fact, a careful inspection of the proof of fliL Theorem 2.1] shows that the non- 
trivial potential fo in (|1.7p obstructs the argument under assumption (|1.3|) . In the 
context of (AO) (i.e. with fo = 0), we can prove that condition (| 1 .8[) is verified 
under assumption (|1.3p . for all Ao <G (0, Aoo). This result (Proposition I3.5[) and 
the spectral conditions (Proposition I3.3[) yield the stability of the standing waves 
i>\{t,x) = c lM u(X){x) 7 for all A G (0^) (Theorem To the best of our 

knowledge, Theorem 13.81 is the first rigorous stability result for the asymptotically 
linear Schrodinger equation. 

Lastly, Section [4] is devoted to an application in nonlinear optics. Following 
our study in of self-focusing planar waveguides in Kerr media, the results of 
Sections [2] and [3] under assumption ( AL) now allow us to discuss the existence and 
stability of TE travelling waves in materials having a saturable dielectric response. 
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1.1. Open problems. Various problems remain unsolved in the higher dimen- 
sional setting. For instance, the global analytic continuation carried out in Subsec- 
tion [52] makes use of the uniqueness of positive even solutions to ()SNLS[) (for each 
fixed A). As far as we know, the uniqueness problem in higher dimension is still 
open. Moreover, due to a lack of compactness coming from the unboundedness of 
the domain, the analytic continuation seems hard to obtain without uniqueness. 

In dimension TV ^ 3, the orbital stability of standing waves along a local branch of 
solutions — such as that obtained in Theorem l2.3l — follows from [7|, Theorem 1 (b)]. 
Nonetheless, even if a smooth global curve of positive radial solutions existed in the 
radial case, our method to prove stability along the whole curve might not work 
in dimension TV > 1 because of an extra term of the form ^^-u'(r) in the radial 
(higher dimensional) version of (|SNLS[) , coming from the expression of Laplacian in 
polar coordinates. We have previously failed to handle this problem in the simpler 
case of the power- type nonlinearity considered in |(|. 

1.2. The prototype. The function / will be required to satisfy numerous struc- 
tural and technical assumptions in order to establish our results. It may be helpful 
to keep in mind the following typical example. 

Example 1.1. Under appropriate conditions on V : R — > R and a > 0, the function 
/ : 1 x R + -> I defined by 

f( x ,s)~V(x) Y ^ (1.9) 

will satisfy all of our assumptions. We will state these conditions in due course, to 
illustrate the general case. 

Terminology and notation. For brevity, we will often refer to properties of 
solutions (X,u) of (|SNLSp . e.g. positivity, evenness etc., while actually meaning 
that u possesses these properties. 

We will work in both the real and the complex Sobolev spaces i? 1 (]R,C) and 
i/ 1 (R, R), depending on whether we consider (|NLSj) or (jSNLSj) . respectively. When 
no confusion is possible, we will merely write H 1 (M.), and similarly for L 9 (R), 
i? 2 (R), etc. All of these spaces will be regarded as real Banach spaces, endowed 
with their usual inner products and norms. 

The symbol C will denote various positive constants, the exact value of which 
does not play an essential role in the analysis. 

2. A GLOBAL CURVE OF SOLUTIONS 

Our approach in this section will be based on previous results about bifur- 
cation for semilinear equations in W N . Firstly, in Subsection l2.1[ we will show that 
a local smooth branch of solutions of (|SNLS|) bifurcates from the line of trivial solu- 
tions at the point (0, 0) € R x i7 1 (R). This will be based on similar results to those 
of Q , where local bifurcation and stability results are established for standing waves 
of the NLS in dimension N 3. Under appropriate symmetry and monotonicity 
assumptions, we will then show in Subsection 12.21 that a version of the implicit 
function theorem can be applied at any positive even solution of (|SNLS|) . thereby 
ensuring global continuation of the local branch. The asymptotic bifurcation results 
of [| will then allow us to discuss the asymptotic behaviour as A — > Aqq. 
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2.1. Bifurcation of small solutions. In [3] we proved local bifurcation and sta- 
bility results for the NLS in dimension N ^ 3. The nonlinearities we considered in 
01 can be written as perturbations — in a sense that will be made more precise 
below — of the signed-power nonlincarity 

g(x,s):=V(x)\s\ p - 1 Sl p>l, (2.1) 

with V £ C 1 (R'' V ). The main hypotheses about g involve a parameter b £ (0,2). 
Roughly speaking, it is required that V(x) ~ \x\~ h as |x| — > oo and that the 
problem be 'subcritical', in the sense that p < 1 + jf-^ ■ We will formulate the exact 
hypotheses in the one-dimensional setting below, but let us already mention two 
differences from the case where N ^ 3. Firstly, if TV" — 1, we must impose b £ (0, 1). 
This is a requirement of the variational formulation of a limit problem involving 
the coefficient \x\~ b . Secondly, in dimensions N = 1,2, the problem is 'subcritical' 
for all p > 1, so we can dispose of the above upper bound. However, we will only 
be interested here in the case of bifurcation from the line of trivial solutions, while 
more general situations are considered in [7j, allowing for asymptotic bifurcation 
(i.e. bifurcation from (0,oo) in R x H 1 (M. N )) as well. This restriction will impose 
another upper bound onp > 1, namely p < 5—26. As can be seen from Example l3.4[ 
this condition is also essential to the stability of the standing waves of (|NLS|K 

We will now state the one-dimensional version of the bifurcation result of ■ It 
is convenient to define 

f(x,s) :=f(x,s 2 ) S , (2.2) 

We then suppose that 

f(x, s) = g(x, s) + r(x, s), (2.3) 
where g is defined in (|2.1[) , V £ C 1 (R) satisfies 

lim \x\ b V(x) = 1 and lim \x\ b xV'(x) = -b, (2.4) 

\x\— >oo >oo 

for some b £ (0,1), and the rest r, defined by (|2.3j) . is 'small' in a precise, technical 
sense. As in [7[, we are dealing here with situations where the linearization of 
(|SNLS|) at u — has purely continuous spectrum. The method we used in Q 
to get bifurcation from the continuous spectrum is by perturbation of the model 
nonlinearity g that was considered earlier in Q. 

A fairly technical method was developped in 0, based on a rescaling and a 
perturbativc argument, using a limit equation involving the nonlinearity g with 
V{x) = |:r|~ b . Thus, via continuation from this limit problem, the asymptotic 
behaviour of V as — > oo turns out to govern the local bifurcation from A = 0. 
Hypotheses about the rest r are formulated in — see 0, (rl)-(r5)] — , ensuring 
that the perturbed nonlinearity retains the main properties of g for small |s|, and 
the same asymptotic behaviour under scaling, in the limit A — > 0. 

In the present context, having (|2.1[) to (|2.4I) in mind, we will formulate these 
assumptions in the one-dimensional setting directly in terms of /. Note that the 
bifurcation analysis for the model nonlinearity g was carried out in 0] in dimension 
N = 1, similarly to the higher dimensional problem treated in Q. 

Let us finally remark that singularities at x = were allowed in . We will 
not need to handle singularities here, and so the present hypotheses are formulated 
in a slightly different manner. 
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(Al) fe C^M 2 ^, /(x, •) G C 2 (M\ {0}) and dj(x, •) G C 1 W for all x £ 
(A2) There exists b G (0, 1), p G (1,5- 26) and s £ (0, ^] such that, 
for < \s\ sC 2s : 

\sd% 2 f{x,s)\ < 

and 

K/tviKjji^'^ ::::: (2.5) 

(A3) Setting 9 := (2 - 6)/(p - 1), we have 

=p|o;|- 6 |s| p - 1 , for alia: 6K\{0}, se 



CM*- 1 , 


M 




1. 


C\x\- h \s\P-\ 


M 




1; 


C\s\P~\ 


\x\ 




1, 


C\x\-»\s\p-\ 


\x\ 




1. 



lim k~ 2; 



lim k- (2+e *> ( - ) dif ( -,k e s) = -b\x\~ b \s\ p - 1 s 1 for all x el\{0}, s e R. 

As usual, we have denoted by d\f and d 2 f the partial derivatives of / with 
respect to its first and second arguments, respectively. 

For brevity, we will refer from now on to the collection of assumptions (AO), 
(Al), (A2) and (A3) as (A). 

Remark 2.1. Note that the conditions f(x, 0) = dif(x, 0) = d 2 f(x, 0) = 0, x G R, 
related to the corresponding conditions on r in 0, (rl)], arc automatically satisfied 
provided that / satisfies (AO). 

Example 2.2. The nonlinearity corresponding to the function / defined by (|1.9j) 
is 

f(x,s) = V(x) 1 ^ °* vel. 

We suppose that V G C 1 (M) and that there exists & G (0, 1) such that V satisfies 
(|2.4[) . and we let a = (p — l)/2. Then assumption (A) is satisfied if < a < 2 — b. 

Theorem 2.3. Let f be defined by (|2.2p and suppose that (A) holds. There exist 
Xq > and a function uq G C 1 ((0, Ao), H 1 (M.)) such that, for all X G (0,Ao), 
(A, uo(A)) is a solution of (jSNLSj) with u (X) G C 2 (R), uo(A) > on R, and 
wo ( A) (x), Mo (A)' (x) — > exponentially as \x\ — > oo. Furthermore, 

lim ||u (A)|| L 2 (R) = lim |ju (A)|| L 2 (R) = lim ||u (A)|| L oo (R) = 0. (2.6) 

A— s-0 A— s-0 A— s-0 

Proof. This follows in the same way as Theorem 1 (a) of Q , using the bifurcation 
results for N = 1 with the model nonlinearity g in Q, rather than those of @ 
dealing with N ^ 3. Note that the condition p < 5 — 2b in (A2) ensures that 
bifurcation occurs from the line of trivial solutions — this is analogous to the 
condition p < 1 + in Theorem 1 of ■ □ 

Remark 2.4. Note that no symmetry or sign assumptions on the nonlinearity are 
required for Theorem 12.31 In fact, the solutions inherit their positivity from the 
sign properties of the limit nonlinearity (|2.1j) . and the local analysis as A — > 0. 



lr This follows from the assumption made in the introduction that / € C7 1 (M X R+,M), but we 
state it here for completeness. 
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2.2. Global continuation. We will now prove that, under appropriate assump- 
tions, the local branch of solutions of (jSNLSj) given by Theorem (23] can be extended 
to a global C 1 curve. In particular, we will now suppose that the problem is sym- 
metric with respect to x = — which will allow us to restrict the discussion to 
the half- line, x G (0,oo) — and that the nonlinearity satisfies some monotonicity 
conditions. Our precise hypotheses are the following: 

(A4) f(-x, s) = f(x, s) for all (x, s) G R x R+; 
(A5) dif(x, s) < and d 2 f{x, s) > for all x, s > 0; 

(A6) (i) <9i/(-,s) G L°°(R) for all s ^ and ||<9i/(-, s)|U=°( R ) is uniformly 
bounded for s in compact subsets of R+ ; 

(ii) <_:/(•, s) £ L°°(R) for all s ^ and {c_/(a:, -)}x6R is cquicontinuous. 
(A7) (AL) holds with /^(O) > Urn /„,(_). 

\x\— >oo 

For brevity, we will refer from now on to the assumptions (A), (A4) to (A7) as 
assumption (A'). 

Example 2.5. In addition to the hypotheses made in Examplc l2.2[ we take a 1, 
and we suppose that V is even, V > on R, and V'(x) < for x > 0. Then the 
function / defined by (|1.9|) satisfies assumption (A'). 

Let us now collect some important consequences of (A'). 

Remark 2.6. 

(a) From (AO), (A5) and (A7), there exists M > such that 

s$ f(x, s) ^ /^(x) ^ M for all (x, s) G R 2 + . (2.7) 

(b) From (AO), (A4) and (A5), /oo is even and non-increasing on [0, oo), with 

lim /oo(s) = 0. (2.8) 

|x|— >oo 

(c) (A6)(ii) implies that ||(?2/(-, s)||ioo( H ) is uniformly bounded and \d?f{x, -)} x ^m. is 
uniformly cquicontinuous on the compact subsets of — see e.g. [171 . Lemma 5.1]. 
Furthermore, it follows by integration that {f(x, -)} x eR is also uniformly equicon- 
tinuous on the compact subsets of R+ . 

(d) The asymptotic linearization (|1.5|) has a principal eigenvalue. Indeed, setting 

, . , J R (u') 2 - foc{x)u 2 dx 
— Aoo '■— im f ^— : , 

Meff^Kjxfo} J R u^ax 

(A7) and (|2.8[) imply Aoo G (0, 00). Furthermore, it follows from the spectral theory 
of Schrodinger operators (see e.g. j23|) that Aoo is the supremum of the spectrum 
of p.5|) . Since hmi3.1_s.oo foa{x) = 0, we have cr oss = (—00, 0], where a ess denotes the 
essential spectrum of (|1.5[) . Hence, Aoo > is the principal eigenvalue of (|1.5[) . 

In order to discuss global continuation, it is convenient to introduce the function 
F : R x i_ x (R) -> fl- 1 (R) defined by 

F(A, u){x) := u"{x) + f{x, u(xf)u{x) - Xu(x), (2.9) 

where i_ _1 (R) denotes the topological dual of H 1 (R), and the right-hand side of 
(|2.12p is interpreted as an element of i_~ x (R) via the canonical identifications: 

(tp, v) H -i xH i = / tpvdx for all (p G L 2 (R), v G if x (R); (2.10) 
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(u",v) H -x xH i =- u'v'dx for all u,v G ff^R). (2.11) 

JR 

It follows easily from assumption (A') that F e C 1 (M.x H^R), H' 1 ^)). We will 
obtain positive even solutions of (|SNLS[) by solving the problem 

F(X,u) = 0, (X,u) e (0,oo) X iT^R), u > on R. (2.12) 

The following result shows that all solutions of (|2.12[) inherit symmetry and mono- 
tonicity from the nonlinearity. 

Lemma 2.7. Let f satisfy (Al), (A4) and (A5), and (A, u) G R x H 1 (R) be a 
solution of (f2~T2]) . Then u(-x) = u(x) for all x ^ 0, u G C 3 (R) with u'(0) = 
and u'(x) < for all x > 0, and (A, it) is a classical solution of (jSNLSj) . 

Proof. First, it is easily seen that weak solutions of (|SNLS[) are in fact classical 
solutions, and it follows from (Al) that they are C 3 . The remaining statements 



then follow by standard arguments — see e.g. the proof of [la, Lemma 2]. □ 



The following lemma establishes further properties of the solutions, in particular 
their exponential decay. We will suppose that (A') holds throughout the rest of 
this section. 

Lemma 2.8. Let (X,u) G R x H 1 (R) be a solution of (f2~T2]) . 

(i) < A < Aoo. 

(ii) For any e G (0, A), let 77 = A — e. Then there exists r e > such that 

\u(x)\ < \\u\\ L „ {m c~^ xl - r '\ for all x G R. (2.13) 

Furthermore, 

l im ^M =± VA. (2.14) 

a;->±oo u(x) 

Proof. To the principal eigenvalue Aoo of (|1.5p corresponds an eigcnfunction (f^ > 
0. The proof of (i) then follows in a similar way to that of 0, Proposition 14 (iv)]. 

Property (ii) was stated in 0, Proposition 14 (hi)] in the context of (|SNLS|) on 
the half-line, but we did not give the proof explicitly there, so we present it here 
for completeness. By (AO), for any e G (0, A), there exists r e > such that 

\x\ ^ r e f(x, s) < e < A Vs^O. 

Define a function 

z(x) := \\u\\ L ^e-^ x \- r '\ xGR, 
and a set il e := {x G R : \x\ ^ r e , z(x) < 0}. For all x G f2 e we have 
u"(x) = [A - f(x, u(x) 2 )]u(x) ^ [A - e]u(x) = nu(x). 

Hence, 

z"{x) = jyiMlioce-^^l-^ - u"(x) 

< »j(||u||i<x,e-v^(l ie l- r «) - u(x)) = r]z(x) Vi G Q e . 

Furthermore, z(x) = ||u||ioo — u{x) 5^ for |x| = r e and li mu i^fy, z(x) = 0. 
Therefore, if fl e ^ 0, it follows by the weak maximum principle jllL Theorem 8.1] 
that z ^ in a contradiction. Hence fl e — and so 

< \\u\\ L °oe-^ x \- r '\ \x\>r e . 
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A similar argument applied to — u yields 

-u(x) ||u||i»e-V5(l x l- r «\ \x\ > r e . 

Since we clearly have < ||u||iooe _ ^l a: l _ ' r «) for |x| < r e , (|2.13p is proved. 

Finally, by de FHospital's rule, 

.. u'(x) 2 2u'(x)u"(x) u"(x) ^2^ x 

hm } ' = lim - ; ' - . v / = hm . . = hm A - f(x,u(x)) = X, 

|a:|->oo U(X) Z |x|-»-oo 2u(x)W (x) |x|^oo U(X) \x\-yoo 

where we have used (AO) and (A6)(ii) in the last equality. Since we know from 
Lemma |2"771 that u is even with v! < on (0, oo), (|2.14l) follows. □ 

Let us now prove that the solutions of (|2.12|) are non-degenerate. This will allow 
us to extend the local curve of solutions obtained in Theorem 12.31 in a smooth 
manner. We will denote by D 2 F the Frechet derivative of F with respect to its 
second argument. 

Lemma 2.9. Let (X,u) elx H 1 (R) be a solution of (|2.12|) . Then the linear 
mapping D 2 F(X,u) : iJ 1 (R) — > i7 _1 (M) is an isomorphism. 

Proof. The linear operator D u F(X,u) : H 1 (M.) —> H^ 1 (M.) is explicitly given by 

D 2 F(X, u)v = v" + [2d 2 f(x, u 2 )u 2 + f{x, u 2 )]v - Xv. 

Note that we can write it as D 2 F(X, u) = R\ + C, where R\v := v" — Xv and 
Cv := [2d 2 f{x,u 2 )u 2 + f{x,u 2 )]v., v G H 1 (R). Since u £ H^R), it follows from 
(AO) and (A6)(ii) that 

lim 2d 2 f(x,u(x) 2 )u(x) 2 + f{x,u{x) 2 ) = 0. (2.15) 

It is then easily seen that C : 7J 1 (R) — > H^ 1 (M.) is a compact linear operator. 
Consequently, since R\v : H 1 (M.) — > i7 _1 (IR) is an isomorphism for all A > 0, we 
need only show that D u F(X,u) is injective. 

To prove this by contradiction, let us suppose that there exists v £ H 1 (M.) \ {0} 
such that 

v" + [2d 2 f(x,u 2 )u 2 + f(x,u 2 )}v = Xv infl-^R). (2.16) 

Clearly, it follows that v € C 2 (R) n H 2 (R) and that the equation holds in the 
classical sense. We will first prove that v is even. Let w £ H 1 (R) be the odd part 
of v, w(x) := \{v{x) — v(—x)), x £ R. We need to show that w = 0. Suppose 
instead that w ^ 0. We have that w £ C 2 (R), w(0) =0 and w is also a solution 
of the linear equation (|2.16p . Without loss of generality, we can then suppose that 
w'(0) > 0. Hence, w > in a right neighbourhood of x = 0. Let xq > be the 
first positive zero of w, or xq = oo if w > on (0, oo). In case Xq < oo we have 
w(xo) = and w'(xo) < 0. 

We now let z := u' . By Lemma \2.7\ z < 0, z £ C 2 (R), and z satisfies 

z" + [2d 2 f(x, u 2 )u 2 + f{x, u 2 )]z + dif{x, u 2 )u = Xz. (2.17) 

In case xq < oo, integrating the Lagrange identity for w and z between x = and 
x = xq, and using w(0) = 0, yields 

z(xo)w' (xq) = / dif(x, u 2 )wudx. 
Jo 
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If X{) = oo, we get 

0= lim z(x)w'(x) — w{x)z'{x) = \ dif(x,u 2 )wudx. 



It follows from the previous discussion that z(xo)w' (xq) > 0. However, in both 
cases the integral on the right-hand side is < by (A5). This contradiction shows 
that we must have w = indeed. Hence, v is even and v'(0) = 0. 

Consequently, integrating the Lagrange identity for u and v over (0, oo) yields 

d 2 f{x,u 2 )u 3 vdx = 0. 

Since d 2 f(x, u(x) 2 )u(x) 3 > for all x > by (A5), it follows that v must have at 
least one zero in (0,oo). Furthermore, for any x > 0, multiplying (|2.16|) by v and 
integrating over [x, oo) yields 

i />oo />oo 

vv" dy + / [2d 2 f(y, u 2 )u 2 + f(y, u 2 )]v 2 dy = A / v 2 dy. 

J X J X 

Hence, integrating by parts, 

poo poo 

v(x)v'(x) =- v'(y) 2 dy+ / [2d 2 f(y, u 2 )u 2 +f{y, u 2 )~\]v 2 dy, for all x > 0. 

J X J X 

But it follows by (|2.15p that there exists r > such that 

2d 2 f(y, u{y) 2 )u 2 + f(y, u(y) 2 ) - A < -A/2 for all y > r, 

and so v'(x)v{x) < for all x ^ r. In particular, there exists x\ > such that 
t;(a;i) = and v{x) ^ for all x > x\. Without loss of generality, we can suppose 
that v'[x\) > and v(x) > for all x > x\. Integrating the Lagrange identity for 
v and z over (xi, oo) yields 

dif(x, u 2 )vu dx. 

In view of (A6), it follows easily from (|2.17l) that z £ H 2 (M.), and so 

lim z(x)v'(x) — v(x)z' ' (x) = 0. 

x— foo 

Therefore, — ^(a;!)?;^^;!) < by (A5), so that v'(xi) < 0. This contradiction finishes 
the proof. □ 

We are now in a position to prove the main result of this section. 

Theorem 2.10. Let assumption (A') hold. There exists u £ C 1 ((0,\ oo ),H 1 (R)) 
such that, for all A £ (0, Aoo), (A,w(A)) is the unique positive solution of (jSNLSI) . 
u(X) £ C 2 (R) n H 2 (R), u(X) is even, and satisfies Lemmas [E7] and fOl 

Furthermore, there is bifurcation from the line of trivial solutions at X = 0, in 
the sense of (|2.6[) . and asymptotic bifurcation at X — A^, in the following sense: 
if A n — > X £ (0, Aoo] as n oo then 

lim \\u(X n )\\ H 2, R) = lim \\u(\ n )\\ L c*>(m = oo <*=^> A = Aoo- (2-18) 

n— >ao n— >oo 

Proof. Global asymptotic bifurcation for (jSNLSp on (0, oo) was established in Q by 
degree theoretic arguments. Corollary 2 of [sj holds under the present hypotheses 
- restricted to the problem on the half-line in the obvious manner. One needs 
only remark that, in hypothesis (f2) of [§], fo can be assumed to be zero, without 
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any change to the proof of Corollary 2. Moreover, under hypothesis (A5), the 
assumption A* := lim sup_j._j.oo foo{x) > fo can be relaxed to an equality, allowing 
for (|2.8j) . By even extension to K, Corollary 2 of 0] then yields a continuous curve 
of solutions of (|2.12[) . (0, Aoo) 9 A — >• u(X), with the asymptotic behaviour (|2.18l) . 

Under hypotheses (A4) and (A5), it can be proved by the method of 'separation 
of graphs', as presented in [32|, that the positive solution of (|SNLS|) is unique, for 
any A £ (0, Aoo). Therefore, the local branch obtained in Theorem 12.31 lies on this 
curve. Thus, the bifurcation behaviour at A = follows from Theorem ,2.31 whereas 
the asymptotic behaviour as A — > Aoo follows from Corollary 2 of [9[. Finally, 
Lemma [2~9l shows that u £ C 1 ((0, Aoo), ^(R)), concluding the proof. □ 

3. Orbital stability of standing waves 

Using the function u £ ^((O, Aoo), H 1 (R 1 K)) given by Theorem [2~10l standing 
wave solutions of (|NLS|) are constructed as 

Mt,x) := e lAt u(A)(x), Ae(0,Aoo). (3.1) 

The mapping A — >• -0a defines a smooth curve of solutions of (jNLSj) in the space 
C([0, oo), ^ 1 (R, C)) nC 1 ((0, co),if- 1 (R, C)). General solutions of i|NLSj) arc func- 
tions ip G C([0, T), H 1 (R, C)) nC 1 ((0, T), ff-^R, C)) satisfying (jNLSj) in the weak 
sense (with the identifications l]2TT0]) - ([2TTI)) ). for all * £ (0,T). Here, T > deter- 
mines the maximal interval of existence of the solution Lp. If T = oo, the solution 
is called global — this is obviously the case for standing waves. 

A prerequisite for the stability analysis of (jNLSj) is the global well-posedness 
of the Cauchy problem. This is thoroughly investigated in for very general 
nonlincarities. We will only need the following result here, which is proved in 
Section 3.5 of 0. 

Theorem 3.1. Let J be defined by (|2~2|) . with f £ C°(M x R+,R) ; and suppose 
that there exist C > and a £ [0, 4) such that 

|/(af, s)[ < C(l + JaDlal, for all (x,s) £R 2 . (3.2) 

For any tpa £ _ff 1 (R., C), there is a unique global solution ip £ C([0, oo), i? 1 (R, C)) 
of (jNLSj) . with initial condition (p(0, •) = fo- 
under assumption (A'), we can take a = in (j3.2|) (see (|2.7j) ). hence the conclu- 
sion of Theorem 13.11 holds . 

Definition 3.2. We say that the standing wave ip\ in (|3.1|) is orbitally stable if 

V e > 3 8 > such that 
for any (global) solution <p(t,x) of (jNLSj) with initial data ipo £ if 1 (R) we have 
\\fo - "(A)Hffi s$ 5 =^ inf 11^(4, •) - e l6 u{X)\\ H i Vi > 0. 

A general theory of orbital stability for infinite-dimensional Hamiltonian systems 
was established in [l2| — see also jig, where this issue was revisited in great detail, 
and applied to (jNLSj) . The stability of a standing wave ip\ is related to spectral 
properties of the linearization of (jNLSj) at ip\. When the spectral conditions are 
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satisfied, it can be inferred from [26| that, in the present context^ 

ip\(t, x) = e lXt u x (x) is orbitally stable if / u(X)(x) 2 dx > 0. (3.3) 

dA J R 

This condition is often referred to as the slope condition. 

3.1. The spectral conditions. Following the discussion in [26| (see in particular 
part (5) of the summary in [26l Section 7.4]), the spectral conditions pertain to the 
linear operators L\,L 2 : H 2 (R,R) C L 2 (K,E) -> L 2 (K,M) defined by 

L\v := -v" - d 2 f(x, u(X))v + Xv, h\v := -v" - /( ^ A)) « + Xv. (3.4) 

Let us denote by M(A) the Morse index of a self-adjoint operator A : D(A) C 
L 2 — s- L 2 , defined by 

M(A) := ^ dim£, 

where £ is the collection of all eigenspaces corresponding to negative eigenvalues of 
A (and M(A) := if £ = 0). We will also denote by o-(A) and a e (A) the spectrum 
of A and the essential spectrum of A, respectively. Then the spectral conditions 
required by the stability analysis are the following: 

(51) inf <r c (L\) > 0, M(L\) = 1 and kcrL^ = {0}; 

(52) inf a e (L 2 ) > 0, = inf a(L 2 ) and kerL 2 = span{u(A)}. 

Proposition 3.3. Let (A') hold, J be defined by ^T^, andu £ C l {{Q,Xoo),H l {R)) 
be given by Theorem\KM Then the operators L\,L\ : H 2 {R) C L 2 (R) -> L 2 (M) 
defined in (|3.4p satisfy (SI) and (S2), respectively. 

Proof. First of all, an argument similar to the proof of [3, Lemma 3.4 (i)] shows 
that all eigenvalues of L\ and L 2 are simple. Also, 

lim d 2 f(x,u(X)) = lim Z^B^l = => inf a e (L\) = X > 0, i = 1, 2. 

We now complete the proof of (SI). First, arguments similar to the proof of 0, 
Lemma 13] show that M(L\) = 1 for A > small enough. Since \\d2f{x, u(X))\\l°° 
depends continuously on A £ (0, Aoo), it follows from the min-max characterization 
of eigenvalues (see e.g. [H, Section XIII. 1]), that the (isolated) eigenvalues of h\ 
depend continuously on A £ (0, Aoo). But kerL^ = {0} for all A £ (0,Aoo) by the 
proof of Lemma 12.91 which prevents eigenvalues from 'crossing zero' as A varies. 
Hence M{L\) = 1 for all A £ (0,Aoo). 

Regarding (S2), noting that 

L\v = -v" - f(x,u(X) 2 )v + Xv, v £ H 2 (R), 

it follows from Theorem 12.101 that u(X) £ kerL 2 . The eigenvalues of L\ being 
simple, kerL 2 = span{w(A)}, A £ (0,Aoo). Finally, since mi a (L%) > 0, is an 
isolated eigenvalue, with corresponding eigenfunction u(X) > 0. It follows that 
= inf er(L 2 ), which completes the proof. □ 



-^Note that the parameter A in [2r3 . Section 7] has an opposite sign from ours. 
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3.2. The slope condition. In order to show that the slope condition is verified, 
we will need the following additional assumption. 

(H) The function 

[x,s) -> — — r 1, X > 0, s > 0, 

d 2 f(x, s)s 

is positive and non-increasing as x increases and s decreases. 

Example 3.4. Under the hypotheses of Examples 12.21 and 12.51 it is easily seen 
that the function / defined by (|1.9[) also satisfies assumption (H), provided that 
—b < xV'(x)/V(x) < for x > 0, and that the mapping x — > xV'(x)/V(x) is 
non-increasing for x > 0. 

Let us now summarize our hypotheses on V : R — > R and a > ensuring that 
the function / defined in Example 11.11 satisfies (A') and (H): 

(a) V £ C^R) is even, V > on R, and V'(x) < for x > 0; 

(b) there exists b £ (0, 1) such that V satisfies (|2.4[) : 

(c) letting p{x) := xV'(x)/V(x), we have that p : (0, oo) — > R is non-increasing, 
with -6 s$ p(x) < for all x > 0; 

(d) Hct<2-6. 

A function V satisfying properties (a)-(c) is given by 

V(x) = {l + x 2 y b / 2 . 

Note that a < 2 — b corresponds to p < 5 — 2b in assumption (A2) (see Example l2.2p . 
Under assumptions (a)-(d), we can take / defined by (|1.9|i in (jNLSp and (|SNLS[) . 
and the conclusions of Theorem 12. 101 and Theorem 13.81 below hold true. 

N.B. A more general form of (|1.9j) can be considered, e.g. f(x,s) = V(x)<f>(s), 
with the above assumptions on V, and cf> £ C 1 (R+) satisfying: 

(01) <f>(s) > and (j)'(s) > for s > 0; 

(4>2) letting $(s) := s<j>' (s) / <)>(s) , there exists a := lim s ^ + $(s) £ [1,2 — b), 

$ : (0, oo) — > R is non-increasing, and < <£>(s) ^ a for all s > 0; 
(03) lim 0(s) = 1. 

In particular, conditions (c) and (02) (simply (c) and (d) in the special case (|1.9[) ) 
ensure that (H) is satisfied. 



In view of Proposition 13.31 and [26( , the orbital stability of the standing waves 
4>\ m p.lj) will be established if we prove the following result. 



Proposition 3.5. Let hypotheses (A') and (H) hold, and u £ C 1 ((0, Aoo), H 1 
be given by Theorem \2.1(A Then 

4r [ u(\)(x) 2 dx > 

for all A £ (0, Aoo)- 



(3.5) 



Our proof of Proposition 13.51 is very similar to that of (1.5) in [g, Theorem 1.7], 
where a power nonlinearity is considered. Since the structure of the nonlinearity is 
quite different in the present case, one needs to check carefully that every step of 
the proof carried out in Section 5.2 of Q works under assumption (A'). Therefore, 
we will present the whole argument here. 
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We will suppose that assumption (A') holds for the rest of this section. We 
already know from Theorem 12.101 — more precisely from bifurcation at A = 
that there is some A > for which (|3.5[) holds. Therefore, by continuity, we need 
only show that 

1 u(X)(x) 2 dx = 2— / u(X)(x) 2 dx^0 VAg(0,A oo ). (3.6) 



dA J a dX 

Now, 



where 



d f°° r°° 
— u(X)(x) 2 dx = 2 u(\)(x)Z(\)(x)dx, 
dA J a Jo 

(JlU 

£(A) := —(A) G H (R), for all A e (0, A*,). 
clA 



We will use a fairly involved integral indentity derived from the equations satished 
by u(A) and £(A), namely (|SNLS|) and 

£" + f{x, u 2 )i + 2d 2 f(x, u 2 )u 2 i = X£ + u. (3.7) 

(We will omit the variables A and/or x when no confusion is possible.) The equation 
for £ is easily obtained by differentiation of the identity ^(A, u(A)) = with respect 
to A, where F is defined in (|2.9j) . We will first prove two lemmas establishing the 
required integral identity and some useful properties of £. Proposition 13 . 51 will then 
be proved using these results. 

The strategy of proof applied here was first used in (l6| . where the authors 
considered a similar problem, however with a power-type nonlinearity (although 
their proof allows for more general situations — albeit not the asymptotically linear 
case) , and a non-trivial linear potential. 

Lemma 3.6. For u,£ £ i? x (R) as above, the following identity holds: 

) POO 

[2f(x,u 2 ) + xd 1 f(x,u 2 )-d 2 f(x,u 2 )u 2 ]utdx = 2X u£dx. (3.8) 

Proof. The proof follows that of 0, Lemma 5.2]. Let us first remark that, using the 
properties of u, it follows from dHTTJ) that £ e C 2 (tt) n iJ 2 (M), and £ is even with 
£'(0) = 0. Integrating the Lagrange identity for (jSNLSj) and (|3.7|) then yields 



u 2 dx = 2 d 2 f(x,u 2 )u A £dx. (3.9) 
On the other hand, multiplying (jSNLSp by u and integrating gives 

/>oo 

f(x, u 2 )u 2 - {u') 2 dx = X u 2 dx. (3.10) 



Now multiplying (|SNLS[) by xu' and integrating by parts (using the exponential 
decay of u') yields 



xf(x, u 2 )uu' — (u 1 ) 2 — xu'u" dx = X j xuv! dx. 
Using (jSNLSP and integrating by parts, it follows that 



oo 



2xf{x, u 2 )uu' — {u') 2 dx = 2X / xuv! dx = 2X I x(^w 2 )'da; 



oo 



II 







oo 


poo 


A 


xu 2 













Jo 



u 2 dx 



poo 

. / u 2 dx. (3.11) 
Jo 
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Furthermore, computing 

d / f u2 \ f u2 f u2 

— (x / /(x,s)ds) = / /(x, s) ds + x/(x, u 2 )2uu' + x I d\f{x 1 s)ds, 

v Jo ' JO JO 

we can substitute the first term of the LHS of (13. lip and we get 

2 2 

(x ^ f(x, s) ds) - [f(x, s) + xd x f{x, a)] ds - {u 1 ) 2 ) dx 

poo 

/ u 2 dx. (3.12) 
Jo 



= -A 

But the first term in the LHS of p,12p can be integrated and yields 

x / /(x, s) ds = lim x / /(x, s) ds = 
Jo o x->oo J Q 

since, by (|2 . 7[) and the exponential decay of u, 

pu(x) 2 

x / /(x,s)ds ^ Mi«(i) 2 -> as x — > oo. 
Jo 

Hence we finally have 

/ / [f(x, s) + xd 1 f(x, s)} ds + (u 1 ) 2 dx = X u 2 dx. (3.13) 
Jo Jo Jo 

Now adding ([3~TU]) and (j3~T3) yields 

2 

poo pu poo 

/ f{x,u 2 )u 2 + / [/(x,s) + xdi/(x,s)]dsdx = 2A / u 2 dx. (3.14) 
Jo Jo Jo 

Next, differentiating (|3.14p with respect to A, we obtain 

[2f(x,u 2 ) + xdif(x,u 2 )+d 2 f{x,u 2 )u 2 }u£dx= u 2 dx + 2A / <dx. 



(3.15) 

Finally, using (J31J) to substitute the first term of the RHS of (j3TT5j) , we get $£E$i . □ 

The following lemma establishes useful properties of £. 

Lemma 3.7. For all A e (0, Aoo), the function £ = £(A) has the following proper- 
ties: £(0) > and t/iere exists a unique xq = xo(A) € (0,oo) suc/i i/iai £(xo) = 0, 
£(x) /or aZZ a; € (0,2:0), and £(x) ^ /or a/Z x £ (xo,oo). 

Proof. We fix A e (0, Aoo), and we simply write £(x) for £(A)(x), x € R. The first 
part of the proof — showing that £(0) > — is easily adapted from that of d, 
Lemma 5.3], using the sign of dif(x,u 2 ) given by (A5). 

To prove the existence of an xo € (0, 00) (xo = Xo(A)) with the asserted prop- 
erties, it is convenient to establish first that there is no a > such that £(x) ^ 
for all x a. This is the part of the proof where some amendments to that of 
H Lemma 5.3] need to be mentionned. We will prove this by contradiction, so we 
assume that such an a exists. Integrating the Lagrange identity for (|2.17p and (|3.7p 
from x ^ a to 00 yields 

-z{x) 2 (x) = £(x)z' (x) - £' (x)z{x) = J wz + 9i/(2/,u 2 )<dy < 0, x > a. 
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Hence (£/z)' > on [a, oo) and £/z is increasing on this interval. But £/z ^ on 
[a, oo) and so there exists a number L such that 

lim 44 = L < 0. 

aj-s-oo z(x) 

Since £ £ iJ 2 (K), we have £'(x) — > as x — > oo. Furthermore, z'(x) —> as x — >• oo 
by (jSNLSj) . It then follows from de l'Hospital's rule that 

L=ii m ^4 =lim m =lim m, 

x-Hx> Z[x) x-yoo Z (x) z-s-oo z"(x) 

as long as the last limit exists. Now 

C \£ + u- f(x,u 2 )^-2d 2 f(x,u 2 )u 2 ^ 



z" Xz — /(x, u 2 )z — <9i/(x, u 2 )u — 2d2f(x, u 2 )u 2 z 
_ Ag/z + u/z- f(x, u 2 )£,/z - 2d 2 /(x, u 2 )u 2 ijz 
A — /(x, u 2 ) — dif(x, u 2 )u/z — 2d%f(x, u 2 )u 2 ' 

where, as x — > oo: 

i/z^-L, u/z^-X' 1 ' 2 , 
/(x, u 2 ) -> by (AO), d 2 f(x, u 2 )u 2 -> by (A6)(ii), 

and 

9i/(x,u 2 ) -> by (USD . 

Hence 

K->oo Z (X) A 

a contradiction. Therefore, there is no a > such that £ (x) ^ for all x ^ a. The 
remainder of the proof is easily adapted from that of [8|, Lemma 5.3], using again 
9i/(x, u 2 ) < 0. We leave the details to the reader. □ 

Proof of Proposition \3.5\ We prove Q3.6P by contradiction. Hence we suppose that 
J °° u(A)(x)£(A)(x) dx = for some A <G (0, Aoo). We omit the dependence on A 
from u and £ for the rest of the proof. It follows from Q3.8P and our assumption 
that 

' 2/(x, u 2 ) + xd lf (x, u 2 )_ 1 u2)u ^ ^ = Q 



S 2 /(x,U 2 )u 2 

Letting 

2/(x ,^ 2 )+x9 1 /( x,u 2 ) 
(9 2 /(x, u 2 )u 2 

it follows from (H) that £ is positive and non-increasing on (0,oo). Using the zero 
xq of £ given by Lemma 13. 71 we can rewrite the identity (|3.16p as 



C(x) := JK iZ . 2\ 2 ~ 1 ' *>°> 



[C(x) - C(x )]a 2 /(x, u > 3 £ dx + C(x ) / 8 2 f(x, u 2 )u 3 d dx = 0. 
By (|3.9p . this becomes 



[C(x)-C(x )]d 2 /(x, W V£dx+^ j[ u 2 dx = 0. 

Since J °° u 2 dx > 0, the properties of £ yield the desired contradiction. The propo- 
sition is proved. □ 
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Theorem 3.8. Suppose that hypotheses (A') and (H) hold. The standing waves of 
(jNLSj) defined by (|3.1J are orbitally stable. 



Proof. In view of the discussion of orbital stability for (|NLS[) in [26[ (see part 
(5) of the summary in Section 7.4 of [2(i]). the result follows immediately from 
Propositions 13.31 and 13.51 □ 



4. Self-focusing planar waveguides 

In this last section we will briefly present an important application of our results 
to nonlinear optics. Self-focusing planar waveguides have been thoroughly inves- 
tigated, both from the physical and the mathematical standpoints. Amongst the 
wide literature about nonlinear waveguides, let us mention 0, 0, Gl El, ED, Hi 



regarding the physics, and [8j, [20|, |22|, |25|, [2JJ, |29| for some mathematical results. 



Historically, the mathematical study of (|NLS|) has grown in parallel to — and was 
largely motivated by — the development of nonlinear waveguide theory. A striking 



example of this close interaction is the early paper [331 ] . where the slope condition 
(|3.3[) was first introduced to discuss the stability of travelling waves in a cylindrical 
waveguide with a saturable dielectric response. (We will precisely be interested in 
saturable materials here, see (|4.2p below.) 

The mathematical modelling of a planar self-focusing waveguide is summarized 
in [§] , where we used bifurcation and stability results similar to Theorem 12.101 
and Theorem 13.81 to discuss the behaviour of TE travelling waves in a waveguide 
with a power-type dielectric response. The planar waveguide is idealized as a slab 
of dielectric material parallel to the xz-plane, having infinite extension in the a;- 
direction, and semi-infinite extension along z, in the z > direction. We will look 
for electromagnetic waves in the optical regime, travelling along the z > half-axis, 
and so x will be transverse to the direction of propagation. According to Maxwell's 
equations, the propagation of a light beam depends on the material, characterized 
by a dielectric response e > (or, alternatively, the refractive index n = y/s). In a 
nonlinear medium, the dielectric response can be decomposed as 

s(x, s) = e L (x) + e NL {x 7 s), 

where Sl{x) and Enl(x, s) respectively denote the linear and nonlinear contribu- 
tions to the dielectric response. As we shall see below, the variable s is proportional 
to the squared modulus of the electric field of the light beam. The dependence 
on the variable x accounts for a medium which is inhomogeneous in the trans- 
verse direction. We will consider the case where the linear contribution, defined by 
El = £(x, 0), is a positive constant, i.e. the material is homogeneous when the beam 
is switched off. On the other hand, we will suppose that Snl(—x, s) = Snl{x, s), 
and snl decreases away from x — 0. This behaviour helps to focus the waves 
around x ~ since, according to Snell's law, the light beam bends towards regions 
with a higher refractive index. A dielectric medium is called self-focusing when 
£nl(x,s) is an increasing function of the variable s (hence of the beam's intensity). 
In [8J] we treated the case of a Kerr medium, where 

e NL {x, s) = a{x)s (4.1) 

for some positive even function a G C 1 (M), decreasing for x > 0. Assumption (|4.1[) . 
leading to a cubic nonlinearity, gives a good approximation, in the low power regime, 
of so-called 'Kerr materials'. However, other materials — e.g. photorefractive 
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materials (see — present a saturation phenomenon as the power of the beam 
becomes large. Namely, 

£nl(x, s) — > e ao (x) as s — > oo, (4-2) 

where s oa (x) adds up to e^, yielding the asymptotic dielectric response + e 00 (x), 
in the limit of high power beams. As can be seen from the equations below, this 
leads to an asymptotically linear nonlinearity (sometimes referred to as a 'saturable 
nonlincarity'). We are now able to deal with this case. 

A TE travelling wave is a special solution of Maxwell's equations in the waveg- 
uide, having an electric field transverse to the sz-plane, of the form 

E{x, y, z, t) = Re (0, U{x)S kz ~ ut) , 0) . (4.3) 

The cnveloppc U of the electric field is such that U e i? 1 (K,IR) so as to verify 
the 'guidance conditions', namely decay of the electromagnetic field at infinity and 
finiteness of the energy density — , and satisfies the Helmholtz equation 

U" + {±) 2 [e L +e NL {x,\U 2 )]U = k 2 U, x e R. (4.4) 

Since we only consider monochromatic waves, i.e. with u > a fixed frequency in 
the optical regime, the time variable t in (|4.3I) does not play an essential role in 
the analysis. The stability of the TE travelling waves (|4.3j) is usually discussed in 
the context of the paraxial approximation, with respect to the larger class of TE 
modes, which are more general solutions of Maxwell's equations than (|4.3p (still 
with an electric field transverse to the a;z-plane). After a rescaling of the variables, 
this approximation leads to the following equation, governing the behaviour of TE 
modes in the waveguide: 

id^ + dlj> + (±) 2 e NL (x, \\i>\ 2 )^ = 0, z > 0, x e R. (4.5) 

A detailed discussion of the paraxial approximation is given in where we also 
explain the reduction of the problem to the canonical form (|4.5p , from which the 
constant ej, has disappeared. In this approach, the TE travelling waves (|4.3[) can 
be approximated by standing wave solutions of (|4.5[) . Since (|4.5p has the form 
of (|NLSp . and (|4.4p that of (jSNLSp . we can apply the results of the previous 
sections, yielding existence of solutions for (|4.4p . hence of TE travelling waves, 
and orbital stability of standing waves of (|4.5p . As explained in 0, Section 6.3], 
a global bifurcation result such as Theorem 12.101 enables one to make the paraxial 
approximation as accurate as desired. Hence stability of the travelling waves (|4.3j) 
can be inferred from Theorem 13.81 — the physical meaning of 'stability' for the 
TE travelling waves (|4.3j) is discussed in Section 6.2 of [lH, where Kerr media arc 
studied using the results in [lj| [ill [24[ ■ 

The rigorous proof of stability of the travelling waves (|4.3p . thus obtained from 
the analysis of (|NLSp . puts on firm mathematical grounds the phenomenon of self- 
trapping, known to physicists since the early 70 's (first predicted theoretically and 



later observed experimentally — see 21|). This phenomenon is briefly described 



as follows. In a self-focusing medium, the nonlinear contribution £jvx(a:, s) is an 
increasing function of s > 0, typically larger near to x = 0. By increasing locally, 
around x = 0, the dielectric response, the light beam induces its own waveguide, 
forcing the waves to focus in the x = region, along the z-axis. When this focusing 
effect balances the dispersion due to the term d 2 x ip in (|4.5p . self-trapping occurs, 
yielding stable guided waves. 
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Our main results for self-focusing planar waveguides with a saturable dielectric 
response (|4.2j) follow from a similar discussion to that in [H, Section 6.3] for the 
Kerr medium (|4.1[) . They can be summarized as follows. Letting 

f(x,s) = (^) 2 e NL (x,ls), 

we suppose that assumptions (A') and (H) arc satisfied. Then, applying Theo- 
rem [5TTU] and Theorem 13. 81 to (|4.4[) and (|4.5[) . we get the following results. 

(1) There exist guided TE travelling waves (|4.3p . corresponding to solutions 
U = Uk of (|4.4p . for all wave numbers fc <E (fci, £3), where 

fci — ^V^ and fc 3 := ((7) 2 £l + Aoo) 1/2 - 

(2) The TE travelling waves P~3)) are stable, for all k g (fci, fc 3 ). 

(3) Defining the power of the beam associated with Uk as 

'=P- [ U k (x) 2 dx, 
ZUJ Jm. 

it follows from our bifurcation results in Theorem 12.101 that 

lim P(fc) = and lim P(k) = 00. 

k — yhi k — yk$ 

Hence, we obtain guided waves from arbitrary low to arbitrary high power. 
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